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ANALYTICAL STUDY OF EFFECTS OF 


SEVERE TURBULENCE ON FLIGHT MOTIONS OF A 
TYPICAL SUBSONIC JET-TRANSPORT AIRPLANE 

By William D, Grantham and Mary S. Adams 
Langley Research Center 

SUMMARY 

An analytical investigation has been made to determine whether the inherent rigid- 
body, uncontrolled stability and response characteristics of subsonic jet-transport air- 
planes could result in a gross upset when the airplane was flying in severe turbulence; 
and, if so, what techniques might aid in the prevention of possible upsets from this cause. 
The investigation consisted mainly of calculations of the motions of an airplane when 
flown through a number of different samples of severe turbulence. 

The results indicated that when the airplane was trimmed at the cruise condition 
and flown through severe turbulence with no pilot inputs, the airplane could become 
grossly upset. The cause of the gross upset encountered in this study was the loss in 
effective dihedral with increase in speed above that for high-speed cruise. The apparent 
sequence of events was that turbulence would sometimes cause a minor upset, generally 
lateral, which caused the Mach number to increase to the point that the effective dihedral 
decreased sufficiently so that the airplane would not raise its low wing (this is the case 
when spiral stability becomes neutral or negative). In the banked attitude, with the cock- 
pit controls remaining fixed in the original position, the airplane would not have enough 
vertical component of lift to support its weight and would therefore continue to lose alti- 
tude and increase airspeed. Several procedures, any one of which could prevent the air- 
plane from being upset in this manner, were identified as follows: 

(a) Decrease altitude and airspeed for turbulence penetration. 

(b) Keep the wings level with loose attitude roll control or with the roll autopilot. 

(c) Use a yaw damper or roll damper. 

The effects of the factors under (b) and (c) are not cumulative, but the use of more than 
one of these does no harm. Control of pitch attitude alone was not adequate for preventing 
the airplane from becoming upset, but it was better than allowing the airplane to fly hands 
off. Since the airplane did not in any case experience a gross upset when the wings were 
kept relatively level, the uncontrolled rigid-body longitudinal response characteristics 
did not seem to be a cause of gross upset for the particular configuration and turbulence 
time histories studied. 



INTRODUCTION 


In the past few years, a number of accidents or near accidents involving subsonic 
jet-transport airplanes have occurred that appear to be related to encountering severe 
turbulence under instrument flight conditions. Such incidents have occurred with several 
different makes and models of subsonic jet-transport airplanes — both commercial and 
military. Some of these incidents occurred when turbulence was encountered at high alti- 
tudes and relatively high airspeeds, while others have occurred at much lower altitudes 
and airspeeds. Typically, the aforementioned incidents were characterized as follows: 

(1) Severe turbulence was encountered while the airplane was being flown on 
instruments. 

(2) The airplane experienced normal accelerations as high as ±3g units 

(29.5 meters/second2), generated by turbulence or by the piloting technique, or by both. 

(3) The cockpit instrument readings became unreliable as a result of rapid fluctua- 
tions of attitude, airspeed, rate of climb, and altitude. 

(4) Airframe structural modes were excited and interfered with the pilot’s normal 
scan of the instruments. 

(5) The pilot's orientation was disrupted. 

(6) The airplane entered a steep dive. 

Practically all incidents reported were encountered during instrument flight conditions, 
and only in those cases where the airplane came out of the overcast and the pilots were 
able to orient themselves visually were recoveries effected and crashes averted. 

Various studies of the turbulence problem have been conducted in the past. In addi- 
tion to the usual investigative studies of the various incidents experienced during turbulent 
flight, other investigations have included computer simulation studies, human centrifuge 
tests, and flight tests. For example, see references 1 and 2. These previous studies, as 
well as the analyses of flight in turbulence presented in reference 3, have led to a better 
understanding of the problems of operating jet-transport airplanes in areas of severe 
turbulence and have particularly highlighted such factors as attitude indicator readability 
and interpretation, turbulence penetration speeds, and control problems such as high stick 
forces and the lack of elevator control power associated with the stalling of the horizontal 
stabilizer drive. 

The present investigation was made to determine whether the inherent rigid-body, 
uncontrolled stability and response characteristics of subsonic jet-transport airplanes 
could result in a gross upset when the airplane was flying in severe turbulence; and, if so, 
what techniques might aid in the prevention of possible upsets from this cause. The 
investigation consisted mainly of calculations of the motions of an airplane when flown 
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through a number of different samples of severe turbulence. The aerodynamic input data 
used in the calculations were based on wind-tunnel tests of a model which was representa- 
tive in general configuration of the class of swept-wing subsonic jet transports. The 
analysis was made by digital computing techniques and included the nonlinearities of 
aerodynamic characteristics and atmospheric properties. 

SYMBOLS 

The units for the physical quantities used herein are presented in both the U.S. 
Customary System of Units and the International System of Units. 

a n normal acceleration, g units (meters/second^) 

ay lateral acceleration, g units (meters/ second^) 

b wing span, feet (meters) 

c mean aerodynamic chord, feet (meters) 

Cj /2 cycles required to damp to one-half amplitude 

rolling-moment coefficient 

C m pitching- moment coefficient 

C n yawing-moment coefficient 

Cx longitudinal-force coefficient 

Cy side-force coefficient 

C% vertical-force coefficient 

g acceleration due to gravity, feet/ second^ (meters/ second2) 

h altitude, feet (kilometers) 

moments of inertia about X, Y, and Z body axes, respectively, 
slug-feet^ (kilograms-meters2) 
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speed of sound at computed altitude, feet/second (meters/second) 

Mach number 

maximum operational Mach number 
mass of airplane, slugs (kilograms) 

rolling, pitching, and yawing angular velocities, respectively, radians/second 

period, seconds 

wing area, feet2 (meters2) 

roll time constant, seconds 

time to damp to one-half amplitude, seconds 

thrust, pounds (newtons) 

components of airplane velocity with respect to inertial space projected along 
X, Y, and Z body axes, respectively, feet/second (meters/second) 

components of airplane resultant velocity Va along X, Y, and Z body 
axes, respectively, feet/second (meters/second) 

incremental linear velocity components along X, Y, and Z body axes, 
respectively, due to turbulence, feet/second (meters/second) 

velocity components of turbulence referenced to X, Y, and Z earth 
axes, respectively, feet/second (meters/second) 

resultant aerodynamic velocity of airplane, feet/second (meters/second) 

maximum operational velocity, knots 

angle of attack, degrees 

angle of sideslip, degrees 
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6 a aileron deflection, positive for right roll command, degrees 

6 e elevator deflection, positive trailing edge down, degrees 

6 S stabilizer deflection, positive trailing edge down, degrees 

6 r rudder deflection, positive trailing edge left, degrees 

£ damping ratio 

9 pitch attitude, degrees 

p air density, slugs/foot3 (kilograms/meter3) 

4> angle of bank, degrees 


\jy angle of yaw, degrees 

o> n undamped natural frequency, radians/second 
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A dot over a symbol represents a derivative with respect to time. 


DESCRIPTION OF AIRPLANE 


The configuration used in this investigation is considered to be representative of 
current subsonic jet-transport airplanes and to be flying at a weight of 175 000 pounds 
(778 435 newtons). The aerodynamic characteristics of the configuration are presented 
in table I. In general, the static aerodynamic derivatives were obtained from refer- 
ence 4, and the dynamic derivatives from unpublished data, with the exception of C m q 
which was obtained from reference 5. The longitudinal- and lateral-directional dynamic- 
stability characteristics of the airplane, without stability augmentation, are presented in 
table n. These data presented in table II are based upon classical linearized, three 
degree-of-freedom equations of motion, and the characteristics presented were derived 
from root solutions of the characteristic equations. The low-speed and high-speed buffet 
boundaries of this airplane, as presented in reference 3, for a gross weight of 
175 000 pounds force (778 435 newtons) and load factors of l.OOg units 
(9.8 meter s/second2) and 1.50g units (14.7 meters/second2), are presented in figure 1. 

PROCEDURES AND CALCULATIONS 


Time histories of the flight motions of the representative subsonic jet-transport 
airplane flying in severe turbulence were calculated by a high-speed digital computer 
which solved the equations of motion and associated formulas listed in the appendix. The 
equations of motion are Euler's equations representing six degrees of freedom along and 
about the airplane body-axes system. (See fig. 2 for illustration of body axes.) 

The initial flight conditions used for the calculations were the efficiency cruise con- 
ditions for this subsonic jet-transport airplane; that is, an altitude of 40 000 feet 
(12.19 kilometers) and a true airspeed of approximately 470 knots (M ~ 0.82). The varia- 
tions of atmospheric properties with altitude were included in the calculations. The 
longitudinal- and lateral-directional static aerodynamic data were used as a function of 
angle of attack and Mach number. The effects of a Mach trim compensator were included 
in the calculations to trim the airplane at the desired Mach number. 

The turbulence used in the present study was measured during severe storm pene- 
trations that were made in conjunction with the National Severe Storms Project (refs. 6 
and 7). The data used consisted of the vertical component w g ' and lateral component 
Vg' of the atmospheric turbulence encountered during four flights through two separate 
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storms. The longitudinal component Ug' was not measured during the storm traverses 
because of instrument failure. However, some longitudinal gusts were arbitrarily used 
in the early stages of the present investigation and it was determined that Ug' had little 
or no effect as to whether the aircraft was upset. (The magnitudes of the longitudinal 
gusts Ug’ were similar to those presented in fig. 3(a) for the lateral gusts Vg\) The 
effects of arbitrary variations in the vertical and lateral gusts were also studied. The 
gusts disturbed the airplane through changes in angles of attack and sideslip. 

The measured turbulence samples used for the computations in this study are pre- 
sented in figure 3 as plots of Wg* and Vg' against time. These gust velocities are 
presented in relation to an earth- axes system, and the equations used to transfer these 
gust components to the airplane body axes are presented in the appendix. (See fig. 2 for 
relationship of airplane body axes to earth axes.) In addition, the gusts were assumed to 
be acting through the center of gravity of the airplane; that is, the airplane was considered 
to be adequately represented by the point mass simulation technique. 

Two control schemes were used briefly during the present investigation in an effort 
to control the motions of the airplane with the aerodynamic control surfaces — a loose- 
attitude or on-off control technique and the more conventional proportional control 
technique. 

The effects of arbitrary variations in some of the more pertinent aerodynamic 
derivatives were also studied. For example, the damping derivatives in roll Qp, pitch 
Cmq> and yaw Cn r , as well as the effective- dihedral parameter and the longitudinal 

static- stability parameter C mQ! , were varied. 

CURRENT RECOMMENDED PILOTING PROCEDURES 
FOR FLYING IN TURBULENCE 

Pilots are aware that when flying in severe turbulence it is possible to impose 
excessive structural loads on the airplane, and that the airplane attitude may reach unde- 
sirable extremes. The flexible sweptwing and high wing loading of current subsonic jet- 
transport airplanes make it probable that any structural damage which might occur in 
severe turbulence will be the result of airplane upset and recovery maneuvers in combina- 
tion with the turbulence, rather than the effects of turbulence alone. Turbulence- 
penetration procedures have therefore been established to minimize airplane attitude 
excursions and thus maintain structural loads within acceptable limits. 

In general, the piloting procedures for severe turbulence penetration are the same 
for all jet-transport airplanes; however, the following procedures are the ones recom- 
mended in reference 3 for the airplane used in the present study: 
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(1) Do not operate above the specified normal cruise ceiling for the airplane weight. 

(2) Maintain indicated airspeed (IAS) of 280 knots or Mach number M of 0.80, 
whichever is lower. Maintain speed by maintaining attitude. 

(3) Maintain control of the airplane with the elevator. After establishing the sta- 
bilizer trim setting for penetration speed, do not change stabilizer trim. 

(4) Make an initial thrust setting for the target speed (IAS = 280 knots or M = 0.80), 
and then do not change thrust unless required by extreme airspeed variation. 

(5) Set the flight director to the desired heading and zero the horizontal bar at the 
desired pitch attitude. 

(6) If the autopilot is not used, use the yaw damper in all cases. 

(7) If the autopilot is used, monitor stabilizer trim, attitude, airspeed, and altitude; 
be alert for an inadvertent autopilot disconnect. Also, if the autopilot is used, the 
altitude-hold mode must be left off for all penetrations of severe turbulence. 

It was possible in the present study to confirm or validate several of these 
recommendations. 


RESULTS AND DISCUSSION 

The calculated results are presented as time histories of pitch attitude angle 0, 
angle of bank (p, angle of yaw xp, normal acceleration a n , lateral acceleration ay, 

Mach number M, altitude h, angle of attack a, and angle of sideslip /3. For the time 
histories of the motions in which various control techniques were used, the magnitudes 
and timing sequences of the control- surface deflections are indicated. Although only 
these more pertinent variables are presented, time histories of all the angles, velocities, 
and accelerations in the equations listed in the appendix were obtained. It should be men- 
tioned that for these calculated time histories the Mach number had an artificial limit of 
0.95 since this was the maximum value for which aerodynamic data were available for 
this configuration. Also, table in presents the range of variation and the root mean 
square of the normal and lateral accelerations experienced during the simulated flights 
discussed. 


Response of Airplane to Severe Turbulence 

The time histories presented in figure 4 represent flights through four individual 
turbulence fields and use the turbulence components presented in figure 3. For each 
simulated flight, the airplane was trimmed at the cruise condition (h = 40 000 feet 
(12.19 kilometers); M = 0.82) and flown through each set of turbulence. The investiga- 
tion consisted of an open loop analysis; that is, no pilot control inputs were used. Exam- 
ination of these time histories shows that the resulting motions in figures 4(a) and 4(b) 
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are appreciably different from the motions indicated in figures 4(c) and 4(d). The motions 
indicated in figures 4(a) and 4(b) would be termed upsets; however, the motions indicated 
in figures 4(c) and 4(d) could not be said to be normal flights. Therefore, the term upset 
must be defined. Reference 8 defines two types of upsets, moderate upsets and gross 
upsets, as follows: "A moderate upset is an excursion to a defined degree beyond any 

normal limit and the following would seem reasonable: 

(a) Too slow — not slower than the stall warning. 

(b) Too fast - not faster than Vj^q + 30 knots or Mj^q + 0.03 Mach number. 

(c) Excessive pitch - not beyond 30° noseup or 20° nosedown. 

(d) Excessive roll - not more than 60° bank. 

A gross upset will then be anything beyond these limits.” From these definitions, the 
motions indicated in figures 4(c) and 4(d) might be termed moderate upsets; whereas, the 
motions indicated in figures 4(a) and 4(b) would definitely be called gross upsets. The 
results of the present study are therefore discussed in relation to gross upsets, and the 
turbulence data presented in figure 3(a) are used unless specifically stated otherwise. 

The effects of variation of airplane stability and control characteristics and gust inputs 
from these conditions are examined. 

Some Factors Affecting Flight Motions in Severe Turbulence 

Effects of gust magnitudes .- The time history presented in figure 4(a) was calcu- 
lated by using no pilot control inputs, no stability augmentation, and no autopilot, and is 
considered the basic flight. Examination of this time history shows that the airplane 
enters a steep spiral dive. The nose pitches down as much as 35° (6 = -35°); the airplane 
rolls to the right as much as 70° (0 = 70°); the heading is changing continuously (0 
approaches 500°); the Mach number reaches its upper limit (M = 0.95); the normal accel- 
eration varies from -1.00 to 4.75g units (a n = -2.00 to 3.75g units from normal lg-unit 
flight); and approximately 29 000 feet (8.84 kilometers) of altitude are lost (h varies 
from 40 000 feet (12.19 kilometers) to approximately 11 000 feet (3.35 kilometers)). 

Several questions arise as to why atmospheric turbulence should cause such results. 
For example, why does the airplane enter a dive and why does the airplane roll to such 
large bank angles ? In order to answer these questions, several additional calculations 
were made. 

The time histories presented in figure 4 indicate that the spectrum content and/ or 
the phasing of the atmospheric turbulence being considered can make a difference as to 
whether the airplane will be upset. In order to obtain some insight as to the peculiarities 
of the turbulence presented in figure 3(a) that caused the airplane to be upset, several 
calculations were made wherein these measured gusts were arbitrarily altered. The time 
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history presented in figure 5 was calculated with all conditions remaining the same as the 
basic flight (fig. 4(a)) except that only the vertical gusts w g ' were used; no lateral gusts 
Vg' were used. The resulting motion indicates that the airplane was not upset and only 
relatively small oscillations in pitch attitude, Mach number, and altitude were experi- 
enced. (The angles of bank and yaw remain zero throughout the time history because no 
lateral-directional disturbances were present.) These results indicate therefore that the 
vertical gusts were not the sole cause of the airplane entering a dive. 

Another storm penetration was calculated wherein the inputs were the same as the 
basic flight except that only the lateral gusts v g ’ were used; w g ' was set equal to zero. 
The resulting motion (fig. 6) shows that again only an oscillatory type of motion was 
experienced; the airplane did not enter a dive and did not achieve unusually large bank 
angles in any one direction. Therefore, the lateral gusts were not the sole cause of the 
airplane achieving bank angles on the order of 50° to 70°. Additional calculations were 
made by using the turbulence presented in figure 3(a), but by arbitrarily reducing w g ' 
and Vg', one at a time, by 75 percent in an attempt to understand further the effects of 
the magnitudes of the vertical and lateral gusts. The resulting time histories, although 
not presented, showed that the airplane was upset in each instance. These results, when 
compared with those of figure 4(d), where there were large gust disturbances but the air- 
plane was not upset, would indicate that the phasing between the vertical and lateral gusts 
may be more important than the magnitude of the atmospheric turbulence insofar as 
grossly upsetting the airplane is concerned. 

Effects of variations in aerodynamic derivatives of airplane .- The question arises 
as to the effects of the inherent stability of the airplane on whether an upset will be 
experienced. From the time histories previously discussed, it appears that the tendency 
of the airplane to roll to unusually large bank angles and to enter a spiral dive could be 
a contributing factor to the jet upset in severe turbulence; that is, no upsets were calcu- 
lated unless large bank angles were achieved. 

Because the spiral mode of the basic airplane used in the present study was neu- 
trally stable (see table II), it would not be expected to cause the airplane to diverge in 
roll after being disturbed. However, an examination of the various time histories, in 
which gross upsets were indicated, shows that the airplane is spirally unstable. This 
characteristic led to a detailed study of the aerodynamic derivatives which normally have 
a large influence on spiral stability. 

Effective dihedral: Previous studies (for example, ref. 9) indicate that as Mach 

number increases to near maximum cruise, the effective dihedral of sweptwing airplanes 
decreases rapidly becomes less negative). This decrease in effective dihedral at 

high subsonic Mach numbers is a fundamental aerodynamic characteristic of aircraft of 
the general configuration under study and is caused by airflow separation on one wing; 
that is, during a sideslip at high speed, the leading wing is less swept than the trailing 
wing, and shock waves will form along the leading wing and cause flow separation and 
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reduce lift on the wing. This shock-wave induced-flow separation on the leading wing 
therefore causes a significant decrease in the effective dihedral. The measured values 
of Cip for the airplane simulated in the present study (ref. 4) are presented in table I 

as a function of angle of attack and Mach number ; at low angles of attack, C ^ changes 
from negative to positive values as the Mach number is increased from 0.80 to 0.90. For 
convenience, C ^ was also plotted against Mach number for the angle of attack at trim 
(a = 0.84° at h = 40 000 feet (12.19 kilometers)) and is presented in figure 7; trim Mach 
number for a= 0.84° and h = 40 000 feet (12.19 kilometers) is 0.82. Note that for the 
angle of attack at trim, the values of C i are positive for Mach numbers above about 

r 

0.86. As C ip is reduced below about -0.003 per degree, the spiral mode becomes 
unstable and this instability worsens as C approaches positive values. If for any 
reason, including turbulence, the Mach number should be increased above approximately 
0.86 and if any lateral disturbance such as turbulence should cause a wing to drop, the 
wing would have to be raised by the pilot or by some artificial means. 

Calculations were made to determine the effects of the magnitude of on the 

tendency of an airplane to become upset in severe turbulence; again, no pilot control 
inputs were used. First, C ip was used as a constant value of -0.0034 per degree, 
which is approximately the measured value for the cruise trim condition (a - 0.84° and 
M = 0.82). The resulting time history is presented in figure 8 and shows that the airplane 
was not upset; 8 never exceeded ±7°, the maximum values of cj) were -14° to 25°, the 
Mach number was not above 0.84 or below 0.78, and a maximum of 500 feet 
(0.15 kilometer) of altitude were lost. Additional calculations were made in which vari- 
ous other constant values of C ^ were used and the resulting time histories indicated 
that as the effective dihedral was decreased, the degree of disturbance of the airplane 
from level flight was progressively worsened while the airplane was flown through the 
turbulence shown in figure 3(a). It is of interest to note the results obtained when 
was held constant at -0.0034 per degree (value used for the calculations presented in 
fig. 8) and the magnitude of the gusts arbitrarily doubled. The resulting time history is 
presented in figure 9, and although large changes in attitude were experienced, this motion 
would not be termed a gross upset. 

These results indicate that the magnitude of can have an appreciable effect 

on whether an airplane will become upset in severe turbulence, and could be the sole 
cause of an upset if the pilot does not keep the wings reasonably level with the lateral 
control. Although in actual operations an obvious part of the piloting task is to maintain 
a wings- level attitude, the effect of marked changes in basic airplane characteristics can 
be an insidious factor. Thus, a rather sudden deterioration of C^, given a Mach num- 
ber increase above M 25 0.82 in turbulence, can seriously degrade an already difficult 
control task, partly because such a change is unexpected. It is possible that very few 
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airline pilots are aware or prepared for such a change in handling qualities because they 
are rarely, if ever, exposed to conditions outside of their normal operating envelope. 


Damping in roll: Variations were also made in the magnitude of the damping-in- 
roll derivative , with all other aerodynamic derivatives being the basic values; that 
is, the same values as were used for the basic calculation presented in figure 4(a). As 
shown in table I, the basic value of Cjp was -0.250 per radian. The additional values 
of this derivative for which time histories were calculated were 2,4, and 10 times this 
basic value = -0.500, -1.000, and -2.500, respectively^. The resulting time histories 

indicate that if the value of C^p were on the order of -1.000, which is four times the 
basic value, the airplane would not be upset. The time history for this condition 
^Cjp = -1.000^ is presented in figure 10. 


Damping in yaw: The damping-in-yaw derivative Cn r was arbitrarily varied in 
the same manner as was C^ p ; that is, 2, 4, and 10 times the basic value, which was 
-0.190 per radian. The resulting time histories showed that when Cn r was as much as 
-0.760 per radian, which is 4 times the basic value, the airplane was not upset. (See 
fig. 11.) Also, as was the case for the variations of C^p, further increases in the 
damping did not improve the resulting motion. 


Other derivatives: Several other aerodynamic derivatives, such as C np , Cj r , 
Cmqj an d C mQ! , were also varied from the basic values, and none of these had any effect 
as to whether the airplane became upset. 


Summary of factors affecting the calculated motions .- Although the factors deter- 
mined during the present open-loop study are not to be taken as the only factors, or even 
the main factors, contributing to the jet upset problem, the following remarks indicate 
several factors regarding the inherent stability and response characteristics that might 
influence the motions of the aircraft in severe turbulence: 


(1) The magnitude of the atmospheric turbulence does not necessarily cause the 
upset; instead the results of this study indicate that the phasing between the vertical and 
lateral gusts is more important. 

(2) The inherent loss of positive effective dihedral as Mach number is increased 
above approximately 0.80 apparently can cause the airplane to enter a spiral dive and thus 
initiate the upset. The apparent sequence of events was that turbulence would sometimes 
cause a minor upset, generally lateral, which caused the Mach number to increase to the 
point that the effective dihedral decreased sufficiently so that the airplane would not raise 
its low wing (this is the condition of neutral, or negative, spiral stability). Since the air- 
plane remained in a banked attitude and the longitudinal controls were still set for the 
original lift coefficient, it did not have enough vertical component of lift to support its 
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weight and continued to lose altitude. The airplane would be, technically, in a spiral dive; 
but the term spiral might give the wrong impressions since it would probably not have 
time to make a full 360° turn during the descent. As the speed increased and air density 
increased, the aircraft would gradually develop more lift and would finally pull out, as 
was the case in figure 4(a), if the bank angle did not continue to increase. 

(3) The results of the present investigation indicate that the lack of sufficient 
amounts of damping in roll and/or yaw can contribute to the jet upset, probably in that 
provision of adequate damping reduced the magnitude of the lateral upset which initiated 
the major loss in altitude. 

(4) It is apparent that the foregoing factors become significant only when there is a 
lack of attention to control, for whatever reason, on the part of the pilot. 

Prevention of Upsets in Severe Turbulence 

Variations in effective aerodynamics .- Obviously, the best way to determine how to 
prevent an upset in severe turbulence is to first determine what factor or factors cause 
the upset, and as discussed in the previous section, three factors have been noted as pos- 
sible contributors to the upset. 

The effect of Mach number on effective dihedral is apparently something that will 
have to be accepted for sweptwing airplanes. However, several possible methods can be 
used to keep C ^ from initiating or contributing to upsets in turbulence. The most 
obvious method is for the pilot not to allow the Mach number to exceed approximately 
0.80. Owing to buffet-boundary considerations, this method is one of the present recom- 
mended piloting procedures for turbulent flight and is discussed subsequently herein. 

The next most obvious method of keeping the variation of C/ with Mach number from 
contributing to the jet upset is for the pilot to combat the loss of effective dihedral by 
flying roll attitude; that is, by keeping the wings relatively level < ±20°) with the 
ailerons. The third alternative is to use a augmenter. This method would require 

the development of a sophisticated sideslip sensing device that could be used on opera- 
tional airplanes in turbulent flight to alleviate the negative dihedral effect at transonic 
Mach numbers. 

Another aid towards preventing jet upsets in severe turbulence is to have an effec- 
tive yaw damper. As discussed previously, if the damping- in-yaw derivative C nr is of 
sufficient magnitude, not less than about -0.800 per radian, the upset is less likely to 
occur. For the particular airplane represented in the present study, the desired value of 
C nr could be obtained by using about 20 percent of the rudder authority. One of the pres- 
ent recommended piloting procedures is for the yaw damper to be on during turbulent 
flight. 
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An effective roll damper Qp which would aid in the prevention of upsets in severe 
turbulence was also determined during the present study. A comparison of the time 
histories presented in figures 11 and 10, wherein a yaw damper and a roll damper were 
used respectively, indicates that the roll damper was, in general, as effective in pre- 
venting the upset as was the yaw damper. In order to achieve the amount of Q n used 
in this case (C^ = -1.000), approximately 100 percent of the aileron authority would be 
required. 

A time history was also obtained wherein the yaw damper and roll damper (C nr and 
Qp increased to four times the basic values) were used simultaneously. The resulting 
motion indicated that when both dampers were used, the motion was similar to that 
obtained when the roll and yaw dampers were used separately. 

Autopilot in turbulence .- Whether the autopilot is used or not, the aircraft and air- 
line industries are in agreement with the fact that the yaw damper should be on for turbu- 
lence penetrations to improve the Dutch roll damping and thus reduce the pilot workload. 
Also, it is generally agreed, as stated previously under Current Recommended Piloting 
Procedures for Flying in Turbulence, that, if the autopilot is used, the altitude- hold mode 
should be off for turbulence penetrations. However, industry apparently has not definitely 
decided whether the autopilot should be used in turbulence penetration. Some factors 
favor having the autopilot off; others favor having it on. Some arguments for and against 
the use of the autopilot in turbulence are presented in reference 3. Reference 3 also 
states that it is neither necessary nor desirable to require that the autopilot be off in 
turbulence in that loss of control which might occur without the autopilot might well be 
prevented if the autopilot were used; also, use of the autopilot leaves the pilot free to 
monitor, which in itself is an important safety factor. 

Although no attempt was made to incorporate an autopilot in the calculations made 
during the present study, as mentioned previously, attempts were made to control the 
attitude of the airplane by using either a loose- attitude control technique or a propor- 
tional control technique. 

Loose-attitude control: Time histories were calculated wherein attempts were 
made to control pitch and bank angles within ±5°, both singularly and in combination, with 
the use of the elevator and aileron, respectively. (Again, the turbulence data of fig. 3(a) 
were used.) These calculations were made with the aforementioned yaw damper 
(C nr = 4 X Basic value) on and off. 

With the yaw damper off, a time history was calculated by using the same inputs as 
were used for the basic calculations (fig. 4(a)) except that an attempt was made to control 
the pitch attitude 6 within ±5° with the elevator being deflected ±5° whenever 6 
exceeded these limits. The resulting motion showed that 6 was held within 5° and -10°; 
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however, the bank angle was as large as 70°, Ai// was greater than 400°, the Mach num- 
ber varied from 0.69 to 0.89, and the airplane lost approximately 10 000 feet (3.05 kilo- 
meters) of altitude. From these results, it might be said that the loose -attitude tech- 
nique of controlling pitch is better than flying the airplane hands off (fig. 4(a)); however, 
this technique would not be said to be an effective system for keeping the airplane from 
becoming upset. Actually, such a system would not be expected to be effective if the air- 
plane is allowed to achieve large bank angles. This case was not repeated with the yaw 
damper on because the results presented in figure 11 (computed for yaw damper on and 
no control inputs) show that the bank and pitch angles experienced never exceeded approxi- 
mately ±15° and ±6°, respectively; therefore, if an effective yaw damper is used, the use 
of a pitch- attitude autopilot would not be required. 

The time history presented in figure 12 shows that when an attempt was made to 
control the bank angle within ±5° (with aileron deflections of ±4° when these limits were 
exceeded), (f> was essentially held within ±10°, 6 was never larger than ±7°, the Mach 

number varied from 0.79 to 0.84, and the altitude did not vary more than approximately 
700 feet (0.21 kilometer) from the trim altitude. Figure 12 also shows that the normal 
and lateral accelerations were approximately ±1.00 g unit and ±0.50 g unit, respectively 
(a n = -0.20 g unit to 1.90 g units and ay = -0.45 g unit to 0.50 g unit); thus, loose- 
attitude control of the bank angle would not be expected to induce excessive structural 
loads on the airplane. When the aforementioned yaw damper was used in conjunction 
with this loose- attitude control of <p, the resulting motion was similar to that obtained 
without the yaw damper. 

The time history presented in figure 13 was computed by using the pitch and roll 
loose-attitude control technique, and the results are very similar to those presented in 
figure 12 where only bank angle was controlled. Figure 13 shows that the elevator was 
required only twice during the simulated flight. Again, when the yaw damper was used 
in conjunction with the pitch and roll attitude control, no apparent effect on the resulting 
motion was obtained. 

The results obtained when the loose- attitude control technique is used are summa- 
rized as follows: 

(1) When the wings were kept relatively level, no gross upsets were experienced. 

(2) Controlling pitch attitude alone was better than flying the airplane hands off but 
was not considered an adequate control procedure. 

(3) When the wings were kept relatively level with the ailerons, the addition of the 
pitch control and/ or the use of the yaw damper showed little improvement in the resulting 
motions. 
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Proportional control: Because all results discussed thus far have shown that no 
gross upsets were obtained when the wings were kept relatively level, the only 
proportional -control scheme investigated was directed toward the control of bank angle 
and was accomplished by simply adding a term C j , in the equations of motion. Effec- 

9 6 a 

tively, the term C is determined from — where the magnitude of 


depends on the gain assigned to 6 a /$. Various time histories were calculated wherein 
6 a/0 was varied, and the results showed that the airplane was controlled very well as 
long as the gain was of sufficient magnitude. Figure 14 presents the motion obtained 
where 6 a /<£ was as low as -0.10, and shows that the gross upset was prevented when 
the airplane was flown through the turbulence presented in figure 3(a). If the 
proportional- control technique can be said to represent a roll-control autopilot, then 
these calculated results would indicate the advantage of using the roll mode of the auto- 
pilot when severe turbulence is encountered. 


Use of proper piloting procedures .- The general piloting procedures for severe 
turbulence penetration have been listed previously herein. These procedures, which were 
based on buff et - boundary considerations, were taken from reference 3, and the results of 
the present study do not contradict any of them. In regard to these recommended proce- 
dures, reference 3 also stated that, "Moderate turbulence will usually require operation 
one cruise level (4000 feet) below the normal cruise ceiling. Protection against greater 
than moderate turbulence will usually require operation two cruise levels (8000 feet) 
below this ceiling." A time history was calculated wherein the inputs were the same as 
the basic case (fig. 4(a)), except that the initial conditions were changed to correspond to 
the aforementioned piloting procedures; that is, the altitude was decreased from 
40 000 feet (12.19 kilometers) to 32 000 feet (9.75 kilometers) (tw-o cruise levels), and 
the true airspeed was decreased from 793 feet per second (241.71 meters per second) 
(which corresponds to M = 0.82 at 40 000 feet (12.19 kilometers)) to 772 feet per 
second (235.31 meters per second) (which corresponds to IAS « 280 knots and M ~ 0.78 
at 32 000 feet (9.75 kilometers)). The resulting motion is presented in figure 15, and as 
can be seen, the airplane was not grossly upset. The pitch attitude was no greater than 
-9° to 5°; the bank angle was no larger than approximately ±25°; the Mach number varied 
from approximately 0.76 to 0.83; and less than 2000 feet (0.61 kilometer) of altitude 
were lost. 


This large effect of airspeed with altitude was unexpected because the dynamic char- 
acteristics of the airplane would not be expected to change considerably for such rela- 
tively small changes in altitude and airspeed. However, because a decrease in altitude 
tends to increase the aerodynamic damping of the airplane, and because the present study 
has shown that no gross upsets were experienced unless large roll angles were achieved, 
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calculations were made to determine the effects of altitude on the dynamic lateral- 
directional stability characteristics of the airplane. The change in the dynamic stability 
between 40 000 feet (12.19 kilometers) and 32 000 feet (9.75 kilometers) was such a 
small amount (table II) that this change would not be the reason the airplane was not 
upset when flying through the same storm (fig. 3(a)) two cruise levels (8000 feet 
(2.44 kilometers)) below the normal cruise altitude. 

The most logical explanation as to why the airplane did not experience an upset 
when flown through turbulence at an altitude of 32 000 feet (9.75 kilometers) (fig. 15) is 
the variation of C ip with Mach number. As stated previously, at 32 000 feet (9.75 kilo- 
meters) the initial Mach number was approximately 0.78, as compared to approximately 
0.82 at an altitude of 40 000 feet (12.19 kilometers). As shown in table I, at low angles 
of attack the effective -dihedral parameter decreases rapidly for Mach numbers greater 
than 0.80. 

Several time histories were calculated wherein the variation of C ^ with Mach 
number was arbitrarily varied to validate the hypothesis that the variation C ip with 
Mach number, rather than increases in aerodynamic damping or variations of dp with 
angle of attack, was the factor that kept the airplane from entering a spiral dive when 
flown at an altitude of 32 000 feet (9.75 kilometers) (fig. 15). The following table sum- 
marizes the general conditions and results: 


Figure 

Initial 

altitude 

Initial 

Mach 

number 

%/ M 

relationship 

Result 

ft 

km 

4(a) 

40 000 

12.19 

0.82 

Basic 1 

Upset 

15 

32 000 

9.75 

.78 

Basic* 

No upset 

17 

40 000 

12.19 

.82 

Curve a2 

Upset 

18 

32 000 

9.75 

.78 

Curve B2 

No upset 

19 

40 000 

12.19 

.82 

Curve c2 

No upset 

20 

32 000 

9.75 

.78 

Curve d2 

Upset 


* C used as function of M and a. 

2 dp used as function of M only, but altered to desired break 
in Cipj M curve. 

Calculations were made to determine whether the variation of with Mach number or 

the variation with angle of attack was the primary cause of the results presented in fig- 
ures 4(a) and 15; and to accomplish this, time histories were calculated wherein dp 
was solely a function of Mach number. The variations of presented in figure 16 as 

curves A and B, which correspond to the measured values at the trim angles of attack 
used for the calculation of the motions indicated in figures 4(a) and 15, respectively, were 
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used as inputs. The resulting time histories are presented in figures 17 and 18, respec- 
tively, and indicate that the variation of C ^ with Mach number was the principal factor 
that brought about the aforementioned results (figs. 4(a) and 15). Actually, examination 
of the time histories of figures 17 and 18 reveals that the resulting motions were similar 
to the original results obtained when was used as a function of angle of attack and 

Mach number (figs. 4(a) and 15, respectively). 

The arbitrary variations of C with Mach number, indicated in figure 16 as 
curves C and D, were used as inputs to validate further the implications of these results. 
Curve C was constructed so that the trim Mach number used for the initial conditions, 
where h = 40 000 feet (12.19 kilometers) and M = 0.82, would fall approximately 

t 

0.02 Mach number before the break in the curve for C^. A time history was calculated 
that corresponds to the basic condition (fig. 4(a)), except that the values of indi- 

cated as curve C were used; the results are presented in figure 19. The resulting motion 
indicated that the upset was not experienced. Similarly, the values of presented in 

figure 16 as curve D were arbitrarily chosen so that the trim Mach number used for the 
initial conditions at the lower altitude (h = 32 000 feet (9.75 kilometers) and M = 0.78) 
would fall approximately 0.02 Mach number after the break in the curve for C^. A time 
history (fig. 20) was calculated that corresponds to the motion presented in figure 15, 
except that the values of indicated as curve D were used. Again, the airplane was 

grossly upset, and a spiral dive was experienced. 

The results of the preceding calculations (figs. 17 to 20) indicate that the variation 
of with Mach number can be a contributing factor in jet upsets when severe turbu- 

lence is encountered and that by following the current recommended piloting procedure of 
decreasing the altitude two cruise levels and flying the airplane at the suggested air- 
speeds, the probability of turbulence increasing the Mach number to such an extent that 
Cip will be reduced to the point of causing spiral divergence will be more remote. This 
recommended piloting procedure should always be followed, if practical, when severe 
turbulence is expected. 


CONCLUSIONS 

An analytical investigation has been made to determine whether the inherent rigid- 
body, uncontrolled stability and response characteristics of subsonic jet-transport air- 
planes could result in a gross upset when the airplane was flying in severe turbulence; 
and, if so, what techniques might aid in the prevention of possible upsets from this cause. 
The investigation consisted mainly of calculations of the motions of an airplane when flown 
through a number of different samples of severe turbulence. The results of this study 
indicate the following conclusions: 
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1. When the airplane was trimmed at the cruise condition and flown through severe 
turbulence with no pilot inputs, the airplane could become grossly upset. 

2. The cause of the gross upset was the spiral instability of the airplane at speeds 
only slightly above the efficiency cruise Mach number. The spiral instability resulted 
from the inherent decrease in effective dihedral with an increase in Mach number. 

3. The apparent sequence of events was that turbulence would sometimes cause a 
minor upset, generally lateral, which caused the Mach number to increase to the point 
that the effective dihedral decreased sufficiently to cause a decided spiral instability 
which in turn resulted in the gross loss of altitude and a further increase in airspeed. 

4. The lack of sufficient amounts of damping in roll and/or yaw can contribute to 
the jet upset. 

5. The phasing between the vertical and lateral gusts may be more important than 
the magnitude of the atmospheric turbulence insofar as grossly upsetting the airplane is 
concerned. 

6. By following the current recommended piloting procedure for turbulence pene- 
tration, that is, by decreasing the altitude two cruise levels and by flying the airplane at 
the suggested airspeeds (Indicated airspeed = 280 knots or Mach number = 0.80), the 
probability of turbulence increasing the Mach number to such an extent that the effective 
dihedral will be reduced to the point of causing spiral divergence will be more remote. 

7. The results obtained when a loose-attitude control technique is used are sum- 
marized as follows: 

(a) When the wings were kept relatively level, no gross upsets were 
experienced. 

(b) Controlling pitch attitude alone was better than flying the airplane 
hands off, but was not considered an adequate control procedure. 

(c) When the wings were kept relatively level with the ailerons, the addi- 
tion of the pitch control and/or the use of the yaw damper showed little improve- 
ment in the resulting motions. 

8. If the proportional-control technique can be said to represent a roll-control 
autopilot, these results would indicate the advantage of using the roll mode of the auto- 
pilot when severe turbulence is encountered. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., July 26, 1968. 
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APPENDIX 


EQUATIONS OF MOTION AND ASSOCIATED FORMULAS 

The equations of motion used in calculating the motions experienced upon pene- 
trating severe turbulence are 

p = qr + + + ° 4 r Sr ) + (V + C ' r r ) 


• Iz “ lx pV a^Sc / . pV aSc^ 

q = I Y pr + 2I Y ( Cm + Cm 5e 6e + Cm 8 s 5 s) + 4Iy ° m q q 


•_ 1 X~ 1 Y^ , pV A ^ 


r = — pq + 


IZ 2I Z 


(Cn^/3 


+ C nc 6p + Cnc 6-rl + 


\ pV A Sb . 

') + 41^ ( Cn p p + Cn r r 


n 5 a a + n 6 r r 


ii = -g sin 8 + vr - wq + 


pV A S 

2m 


( C X + C Xg e 6 e + C X 5g 5 s) + h 


• P^a ® / \ P v A ou / \ 

v = g cos 0 sin <P + wp - ur + (C y ^/3 + Cy^a + + _^_^C Yp p + C Yl .rj 


pV A Sb / 


PV a 2 s . 

w = g cos 0 cos <p + uq - vp + ■ ^ — (C z + C Z g 6 e + 


The following formulas were also used: 


V A = (U; 


a 2 + v a 2 + w a 2) 


1/2 



u a = u - Ug 
v a = v - Vg 
W a = w - Wg 

ug = -Ug' cos xp cos 9 - Vg' sin \p cos 6 + Wg' sin 9 
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APPENDIX 


Vg = -Ug'(cos xp sin <p sin 0 - sin xp cos (p ) 

- Vg'(cos xp cos 4> + sin xp sin <p sin 0) - w g 

Wg = -Ug’(cos xp cos (p sin 0 + sin xp sin <p) 

- Vg'(sin xp cos cp sin 0 - cos xp sin <p) - 


, _i w a 
a = tan 1 — 

Uo 


,3= sin' 1 1* 

V A 


M = 


Ya 

K 


0 = q cos (p - r sin cp 


0 = p + q tan 0 sin cp + r tan 0 cos cp 


'P = 


r cos cp + q sin <p 
cos 0 


h = u sin 0 - v cos 0 sin cp - w cos 0 cos 


a -n = ~ 


g = 32.17 


w - uq + vp - g cos 0 cos cp 

g 

3956.66 'y* 

h 


, 3956.66 + 


5280/ 


'(sin <p cos 0) 


'(cos <p cos 0) 
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I 


TABLE I .- TABULATION OF AERODYNAMICS USED IN COMPUTATIONS 


( a ) Derivatives used as function of a and M 


a 

M = 0.40 

M = 0.60 

M = 0.80 

M = 0.90 

M = 0.95 

M = 0.40 

M = 0.60 

M = 0.80 

M = 0.90 

M = 0.95 

M = 0.40 

M = 0.60 

M = 0.80 

M = 0.90 

M = 0.95 




Cm 





Cx 





c z 



-20 

0.185 

0.160 

0.125 

0.090 

0.182 

- 0.034 

- 0.044 

- 0.059 

- 0.065 

- 0.073 

0.60 

0.60 

0.59 

0.58 

0.67 

-16 

.090 

.072 

.045 

.022 

.095 

-.034 

-.044 

-.061 

-.061 

-.072 

.58 

.55 

.51 

.50 

.59 

-12 

.045 

.025 

.010 

0 

.067 

-.034 

-.044 

-.061 

-.056 

-.071 

.51 

.47 

.44 

.42 

.49 

-8 

.035 

.020 

0 

-.022 

.047 

-.034 

-.043 

-.055 

-.050 

-.071 

.38 

.34 

.34 

.31 

.39 

-4 

0 

-.020 

-.040 

-.055 

.015 

-.032 

-.037 

-.043 

-.043 

-.071 

.13 

.11 

.11 

.10 

.22 

0 

-.060 

-.065 

-.087 

-.120 

-.042 

-.026 

-.025 

-.029 

-.048 

-.073 

-.25 

-.27 

-.31 

-.29 

-.16 

2 

-.085 

-.095 

-.120 

-.140 

-.085 

-.017 

-.017 

-.023 

-.051 

-.077 

-.42 

-.45 

-.52 

-.48 

-.39 

4 

-.105 

-.110 

-.125 

-.150 

-.150 

0 

.001 

-.020 

-.053 

-.077 

-.60 

-.63 

-.70 

-.75 

-.60 

6 

-.110 

-.105 

-.130 

-.158 

-.168 

.024 

.020 

-.024 

-.052 

-.071 

-.75 

-.79 

-.79 

-.75 

-.79 

8 

-.110 

-.110 

-.140 

-.172 

-.188 

.032 

.012 

-.032 

-.052 

-.072 

-.88 

-.90 

-.85 

-.84 

-.93 

10 

-.120 

-.145 

-.170 

-.170 

-.190 

.027 

-.010 

-.040 

-.057 

-.067 

-.99 

-.97 

-.91 

-.92 

- 1.01 

12 

-.135 

-.182 

-.220 

-.175 

-.202 

.009 

-.024 

-.048 

-.060 

-.068 

- 1.08 

- 1.04 

-.96 

-.99 

- 1.07 

14 

-.170 

-.230 

-.250 

-.205 

-.235 

-.008 

-.030 

-.055 

-.064 

-.070 

- 1.16 

- 1.09 

- 1.01 

- 1.07 

- 1.15 

16 

-.230 

-.275 

-.270 

-.250 

-.280 

-.018 

-.036 

-.058 

-.066 

-.073 

- 1.22 

- 1.14 

- 1.08 

- 1.16 

- 1.26 

18 

-.330 

-.335 

-.350 

-.315 

-.345 

-.024 

-.041 

-.061 

-.067 

-.074 

- 1.26 

- 1.18 

- 1.14 

- 1.22 

- 1.32 

20 

-.410 

-.385 

-.372 

-.370 

-.388 

-.028 

-.045 

-.060 

-.068 

-.074 

- 1.29 

- 1.22 

- 1.22 

- 1.29 

- 1.38 




% 





C n/3 





C Y P 



-20 

0.0034 

0.0035 

0.0033 

0.0033 

0.0031 

0.0033 

0.0038 

0.0035 

0.0037 

0.0048 

- 0.0215 

- 0.0215 

- 0.0209 

- 0.0213 

- 0.0221 

-16 

.0033 

.0032 

.0033 

.0024 

.0025 

.0033 

.0038 

.0035 

.0037 

.0048 

-.0217 

-.0219 

-.0213 

-.0215 

-.0226 

-12 

.0024 

.0020 

.0020 

.0007 

.0006 

.0033 

.0038 

.0034 

.0037 

.0048 

-.0219 

-.0223 

-.0216 

-.0217 

-.0231 

-8 

.0002 

0 

-.0003 

-.0016 

-.0020 

.0032 

.0037 

.0034 

.0037 

.0047 

-.0220 

-.0225 

-.0218 

-.0220 

-.0234 

-4 

-.0020 

-.0019 

-.0028 

-.0029 

-.0010 

.0032 

.0036 

.0034 

.0037 

.0046 

-.0215 

-.0221 

-.0220 

-.0221 

-.0235 

0 

-.0031 

-.0033 

-.0044 

.0005 

.0005 

.0029 

.0033 

.0032 

.0037 

.0039 

-.0198 

-.0200 

-.0210 

-.0210 

-.0211 

2 

-.0033 

-.0038 

-.0043 

.0013 

.0003 

.0026 

.0028 

.0030 

.0036 

.0031 

-.0190 

-.0186 

-.0203 

-.0200 

-.0195 

4 

-.0036 

-.0041 

-.0022 

.0008 

-.0009 

.0026 

.0023 

.0033 

.0033 

.0029 

-.0187 

-.0180 

-.0200 

-.0190 

-.0200 

6 

-.0036 

-.0041 

-.0026 

-.0029 

-.0031 

.0025 

.0023 

.0036 

.0031 

.0027 

-.0188 

-.0174 

-.0210 

-.0198 

-.0208 

8 

-.0035 

-.0033 

-.0025 

-.0022 

-.0045 

.0025 

.0031 

.0038 

.0033 

.0022 

-.0190 

-.0190 

-.0210 

-.0210 

-.0213 

10 

-.0030 

-.0023 

-.0021 

-.0015 

-.0063 

.0029 

.0034 

.0039 

.0033 

.0017 

-.0190 

-.0195 

-.0213 

-.0208 

-.0198 

12 

-.0024 

-.0020 

-.0023 

-.0032 

-.0046 

.0033 

.0035 

.0038 

.0031 

.0014 

-.0195 

-.0200 

-.0215 

-.0213 

-.0189 

14 

-.0017 

-.0015 

-.0029 

-.0035 

-.0041 

.0037 

.0037 

.0036 

.0033 

.0014 

-.0201 

-.0202 

-.0213 

-.0215 

-.0180 

16 

-.0011 

-.0012 

-.0028 

-.0029 

-.0038 

.0037 

.0037 

.0037 

.0033 

.0012 

-.0190 

-.0200 

-.0208 

-.0218 

-.0185 

18 

-.0012 

-.0014 

-.0023 

-.0033 

-.0037 

.0031 

.0037 

.0043 

.0032 

.0010 

-.0107 

-.0185 

-.0210 

-.0223 

-.0190 

20 

-.0015 

-.0017 

-.0022 

-.0035 

-.0037 

.0028 

.0036 

.0043 

.0031 

.0007 

-.0186 

-.0182 

-.0212 

-.0227 

-.0194 



Cm q 





C _^__ 








-20 

30 

4 

-20 

-40 

-38 

0.172 

0.182 

0.174 

0.168 

0.169 






-16 

8 

-5 

-3 

-32 

-33 

.162 

.172 

.169 

.162 

.164 






-12 

-7 

-15 

-3 

-25 

-29 

.142 

.152 

.153 

.146 

.147 






-8 

-11 

-21 

-19 

-22 

-27 

.113 

.122 

.125 

.120 

.120 






-4 

-12 

-20 

-19 

-24 

-19 

.073 

.080 

.087 

.083 

.083 






0 

-13 

-16 

-20 

-22 

-18 

.020 

.026 

.032 

.032 

.032 






2 

-14 

-16 

-21 

-20 

-20 

-.013 

-.008 

-.002 

0 

0 






4 

-16 

-19 

-23 

-23 

-28 

-.042 

-.041 

-.034 

-.033 

-.035 






6 

-17 

-18 

-13 

-12 

-18 

-.059 

-.058 

-.053 

-.050 

-.056 






8 

-17 

-15 

-15 

-30 

-13 

-.073 

-.072 

-.070 

-.067 

-.070 






10 

-14 

-10 

-7 

-16 

-38 

-.088 

-.087 

-.082 

-.080 

-.084 






12 

-11 

-5 

-3 

-46 

-8 

-.103 

-.101 

-.096 

-.096 

-.100 






14 

-6 

-2 

-10 

-9 

-20 

-.113 

-.110 

-.105 

-.105 

-.104 






16 

15 

4 

-10 

-43 

-29 

-.121 

-.118 

-.110 

-.110 

-.108 






18 

30 

7 

-50 

-27 

-35 

-.129 

-.125 

-.113 

-.114 

-.110 






20 

25 

5 

-17 

-41 

-37 

-.136 

-.131 

-.116 

-.116 

-.112 
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TABLE L- TABULATION OF AERODYNAMICS USED IN COMPUTATIONS - Concluded 


CY «a 


(b) Derivatives used as constants 


0.00093 

0.00030 

- 0.250 

0.200 

- 0.030 

- 0.0125 

0 

- 0.0020 

- 0.190 

0.0025 

0.0010 

0 

0.0048 

- 0.110 

0.375 

- 0.010 

- 0.004 
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TABLE H.- DYNAMIC- STABILITY CHARACTERISTICS OF AIRPLANE 
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TABLE HI.- RANGE OF VARIATION AND ROOT MEAN SQUARE OF NORMAL AND 
LATERAL ACCELERATIONS EXPERIENCED DURING SIMULATED FLIGHTS 

IN SEVERE TURBULENCE 


Figure 

Maximum normal 
acceleration, 
a n , g units 

Maximum lateral 
acceleration, 
a y , g units 

Root mean square 

a n 

a y 

4(a) 

-1.00 to 4.75 

-0.48 to 0.75 

0.891 

0.148 

4(b) 

0.20 to 3.00 

-0.25 to 0.25 

.723 

.083 

4(c) 

-0.50 to 1.60 

-0.35 to 0.30 

.310 

.106 

4(d) 

0 to 2.00 

-0.30 to 0.40 

.249 

.078 

5 

-0.20 to 1.80 


.293 

0 

6 

0.57 to 1.10 

-0.40 to 0.50 

.065 

.115 

8 

-0.20 to 1.90 

-0.40 to 0.50 

.298 

.115 

9 

-0.90 to 2.50 

-0.90 to 1.10 

.606 

.221 

10 

-0.20 to 1.90 

-0.40 to 0.48 

.296 

.101 

11 

-0.20 to 1.80 

-0.34 to 0.37 

.297 

.071 

12 

-0.20 to 1.90 

-0.45 to 0.50 

.298 

.115 

13 

-0.20 to 1.80 

-0.46 to 0.49 

.291 

.115 

14 

-0.20 to 1.90 

-0.40 to 0.50 

.296 

.117 

15 

-0.40 to 2.40 

-0.50 to 0.50 

.379 

.134 

17 

-0.20 to 4.00 

-0.50 to 0.50 

.779 

.149 

18 

-0.60 to 2.30 

-0.40 to 0.40 

.382 

.126 

19 

-0.20 to 1.90 

-0.50 to 0.50 

.304 

.122 

20 

-0.50 to 3.60 

-0.40 to 0.40 

.722 

.134 
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Indicated airspeed (IAS), knots 


Figure 1.- Indication of low-speed and high-speed buffet boundaries of simulated airplane. Weight = 175 000 pounds force (778 435 newtons). 
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Time, sec 

(a) Basic set of gust samples. 

Figure 3.- Indication of magnitude and phasing of four sets of gust samples used in calculations. 
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(b) Second set of gust samples. 
Figure 3.- Continued. 
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(c) Third set of gust samples. 
Figure 3.- Continued. 
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(c) Response to gust samples of figure 3<c). 
Figure 4.- Continued. 
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igure 5.- Indication of motion obtained when only vertical gusts w g * of figure 3(a) were used. 
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Figure 5.- Concluded. 
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feet 


deg 


p, 

deg 


*e. 

deg 


So 

deg 



Time, sec 


Figure 6.- Concluded. 
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Figure 8.- Concluded. 
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Time, sec 


Figure 9.- Time history of motion obtained when was held constant at -0.0034 per degree and gusts of figure 3(a) were 

arbitrarily doubled in magnitude. 






Figure 9.- Concluded. 
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Time, sec 


Figure 10.- Indication of motion obtained when C[ n was set equal to -1.000 per radian and gusts of figure 3(a) were used. 
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Figure 10.- Concluded. 
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0 




Figure 12.- Indication of motion obtained when bank angle was controlled loosely with ailerons. 
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Figure 12.- Concluded. 
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Time, sec 


Figure 15.- Time history of motion obtained when initial altitude and airspeed were 32 000 feet (9.75 kilometers) and 

M = 0.78 ( respectively. 
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Figure 17.- 
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Indication of motion obtained when curve A of figure 16 was used for C in. Initial conditions: h 

(12.19 kilometers) and M = 0.82. 
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Figure 18.- Indication of motion obtained when curve B of figure 16 was used for C in. initial conditions; h = 32 000 feet 

{9.75 kilometers) and M = 0.78. 
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Figure 19.- Indication of motion obtained when curve C of figure 16 was used for C in. Initial conditions: h 

(12.19 kilometers) and M = 0.82. 
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Figure 19.- Concluded. 
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Figure 20.- Time history of motion obtained when curve D of figure 16 was used for C Ir. Initial conditions: h = 32 000 feet 

(9.75 kilometers) and M = 0.78. 
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Figure 20.- Concluded. 
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